Evolution & Development

The first vertebrate mineralized skeleton was an external bony armor coated with dental structures. The subsequent emergence of a mineralized endoskeleton and of teeth are considered key innovations in the diversification of vertebrates. Although time clouds our understanding of the initial evolution of these mineralized structures, recent re-emergences may shed light on the underlying processes. Loricarioid catfishes are a lineage that, much like the ancestral vertebrates, bear denticle-clad bony armor from head to tail. Loricarioid denticles (LDs) and oral teeth are very similar in superstructure. We show here that other extra-oral dental structures are found as ancestral characters only in lineages that are distantly related to loricarioids such as sharks or coelacanth, indicating that LDs have independently re-emerged in loricarioid catfishes. We investigate whether the similarities between LDs and teeth extend to their developmental and genetic context, and how their development compares to that of other vertebrate integument structures. Our detailed study of the development of LDs, and gene expression analyses through in situ hybridization confirm that all 12 genes from the tooth-forming gene regulatory network (oGRN) are expressed in developing LDs in a similar way as they are expressed in developing teeth. We then compare the developmental, structural, and genetic aspects of LD and teeth with that of other integument appendages such as fish scales, shark dermal denticles, feathers and hairs. We find that LDs share all developmental cues with teeth and, to a lesser extent, with the other vertebrate integument structures. Taken together, our results indicate that denticles have re-emerged on the trunk of loricarioid catfishes through the ectopic co-option of the oGRN rather than the resurrection of an ancestral trunk-specific denticle genetic pathway.

Hemichordates are an important group for investigating the evolution of bilaterian nervous systems. As the closest chordate outgroup with a bilaterally symmetric adult body plan, hemichordates are particularly informative for exploring the origins of chordates. Despite the importance of hemichordate neuroanatomy for testing hypotheses on deuterostome and chordate evolution, adult hemichordate nervous systems have not been comprehensively described using molecular techniques, and classic histological descriptions disagree on basic aspects of nervous system organization. A molecular description of hemichordate nervous system organization is important for both anatomical comparisons across phyla and for attempts to understand how conserved gene regulatory programs for ectodermal patterning relate to morphological evolution in deep time. Here, we describe the basic organization of the adult hemichordate Saccoglossus kowalevskii nervous system using immunofluorescence, in situ hybridization, and transgenic reporters to visualize neurons, neuropil, and key neuronal cell types. Consistent with previous descriptions, we found the S. kowalevskii nervous system consists of a pervasive nerve plexus that is concentrated in the anterior, along with nerve cords on both the dorsal and ventral sides. Neuronal cell types exhibited clear anteroposterior and dorsoventral regionalization in multiple areas of the body. We observed spatially demarcated expression patterns for many genes involved in synthesis or transport of neurotransmitters and neuropeptides but did not observe clear distinctions between putatively centralized and decentralized portions of the nervous system. In the trunk, there is a clear division of cell types between the dorsal and ventral cords suggesting differences in function. The plexus shows regionalized structure and is consistent with the proboscis base as a major site for information processing rather than the dorsal nerve cord. The absence of neural processes crossing the basement membrane into muscle and extensive axonal varicosities suggest that volumetric transmission may play an important role in neural function. These data now facilitate more informed neural comparisons between hemichordates and other groups and contribute to broader debates on the origins and evolution of bilaterian nervous systems.

The potential for variation and the capacity to evolve in response to ecological opportunity are important aspects of an adaptive radiation. Identifying the origin of phenotypic variation, in which natural selection might act upon, is a major goal of evolutionary developmental biology. The New World leaf-nosed bats (phyllostomids) are a textbook example of an adaptive radiation. Their cranial morphology is diverse along relative facial length, which is related to their diets. We previously used geometric morphometrics to reveal peramorphosis, a type of heterochrony, in the cranial evolution among phyllostomid bats. We then demonstrated that the mechanism of peramorphic diversity in phyllostomid rostrum length resulted from altered cellular proliferation. Here, we investigate the progenitors of the face, the cranial neural crest, and a key signaling pathway related to their proliferation and differentiation into mature tissues: the bone morphogenetic protein (BMP). With geometric morphometrics, immunofluorescence, and confocal imaging--in three phyllostomid species and one outgroup bat species--we show the molecular patterns that underlie the adaptive and innovative traits seen in phyllostomid bats. Then, with mouse genetics, we mimic the BMP molecular pattern observed in nectar feeding bats and recapitulate the elongated morphological variation in mice. Surprisingly, we also observe an expansion in the nose-tip of mice, akin to the expanding leaf-nose tissue in phyllostomid bats. These data, combined with the mouse genetics literature on BMP signaling, suggest the BMP developmental pathway plays a central role in shaping the craniofacial variation necessary for adaptation in bats. Further, we speculate that the BMP signaling pathway could underlie other bizarre facial phenotypes in mammals that are derived from frontonasal mesenchyme, such as the proboscis. Overall, this study combines a comparative framework to developmental data, with a genetic approach, to directly investigate the role of development on complex morphology.

Fibroblast growth factor (FGF) signaling is central to chordate development and has been extensively remodeled in tunicates. Recent findings show that appendicularians have massively lost all ancestral chordate Fgf subfamilies except two, the Fgf9/16/20 and Fgf11/12/13/14 subfamilies, which in contrast have undergone a burst of lineage-specific duplications and diversification into novel paralogs, in an evolutionary scenario that we have named "Less, but More". Here, we investigate the downstream effects of the Fgf losses and duplications ion Fgf receptors (FgfRs) and intracellular RTK components in the appendicularian Oikopleura dioica. We show that the single ancestral FgfR gene has expanded into three paralogs (FgfRa-c), which are conserved across cryptic O. dioica species, yet highly divergent from other chordates. Despite strong sequence divergence, structural modeling indicates preservation of canonical FgfR architecture. Expression analyses reveal distinct spatiotemporal patterns: FgfRa and FgfRb are maternally supplied and enriched in mesodermal derivatives, whereas FgfRc is restricted to neural and epithelial tissues. Genome surveys of downstream RTK pathways show conservation of core RAS/MAPK, PLC{gamma}/PKC, and PI3K/AKT cascades, but with losses of classical Ras genes and several adaptors, suggesting a lineage-specific simplification of transduction complexes. Transduction gene expression shifts from broad maternal ubiquity to tissue-specific domains, particularly in brain, notochord, muscle, and gonadal primordia throughout embryonic and larval development. Appendicularians appear as the only non-vertebrate chordate lineage that recapitulate the vertebrate-like FgfR expansion following Fgf ligands diversification. Downstream components, however, evolved more conservatively, tending toward simplification, reinforcing the view that appendicularians generate signaling innovation despite extensive gene loss.

Comparing gene regulatory patterns between seemingly similar phenotypic traits can provide important insights on the molecular mechanisms underlying the evolution of those traits. In this study, we investigate the molecular basis of the formation of a spade-shaped caudal fin, which is a rare phenotype among teleost fish characterized by an elongated medial region of the fin. We examined the expression patterns of candidate fin-shape genes in the spade-shaped caudal fin of the related species Lamprologus tigripictilis, an East African cichlid in the tribe Lamprologini. The candidate gene set consisted of a previously identified gene regulatory network (GRN) associated with the elongation of fin regions in another Lamprologini cichlid species and further genes selected on the basis of co-expression data and transcription factor prediction. Unexpectedly, the anatomical features of elongated fin rays differed and gene expression patterns associated with fin elongation were only weakly conserved between the two related species. We report 20 genes and transcription factors (including angptl5, cd63, csrp1a, cx43, esco2, gbf1 and rbpj), whose expression levels differed between the elongated and the short caudal fin regions of L. tigripictilis, and which are therefore candidates for the regulation of the spade-like fin shape.

Two genes, Distal-less (Dll) and spalt (sal), are known to be involved in establishing nymphalid butterfly wing patterns. They function in several ways: in the differentiation of the eyespots central signaling cells, or foci; in the differentiation of the surrounding black disc; in overall scale melanisation (Dll); and in elaborating marginal patterns, such as parafocal elements. However, little is known about the functions of these genes in the development of wing patterns in other butterfly families. Here, we study the expression and function of Dll and sal in the development of spots and other melanic wing patterns of the Indian cabbage white, Pieris canidia, a pierid butterfly. In P. canidia, both Dll and Sal proteins are expressed in the scale-building cells at the wing tips, in chevron patterns along the pupal wing margins, and in areas of future scale melanisation. Additionally, Sal alone is expressed in the future black spots. CRISPR knockouts of Dll and sal showed that each gene is required for the development of melanic wing pattern elements, and repressing pteridine granule formation, in the areas where they are expressed. We conclude that both genes likely play ancestral roles in organising distal butterfly wing patterns, across pierid and nymphalid butterflies, but are unlikely to be differentiating signalling centers in pierids black spots. The genetic and developmental mechanisms that set up the location of spots and eyespots are likely distinct in each lineage.

Hox genes are transcriptional regulators that elicit cell positional identity along the anterior-posterior region of the body plan across different lineages of Metazoan. Comparison of Hox gene expression across distinct species reveals their evolutionary conservation, however their gains and losses in different lineages can correlate with body plan modifications and morphological novelty. We compare the expression of eleven Hox genes found within Streblospio benedicti, a marine annelid that produces two types of offspring with distinct developmental and morphological features. For these two distinct larval types, we compare Hox gene expression through ontogeny using HCR (hybridization chain reaction) probes for in-situ hybridization and RNA-seq data. We find that Hox gene expression patterning for both types is typically similar at equivalent developmental stages. However, some Hox genes have spatial or temporal differences between the larval types that are associated with morphological and life-history differences. This is the first comparison of developmental divergence in Hox genes expression within a single species and these changes reveal how body plan differences may arise in larval evolution.

BackgroundA hallmark of arthropods is their segmented body, and the so-called Drosophila segmentation gene cascade that controls this process serves as one of the best-studied gene regulatory networks. An important group of segmentation genes is represented by the pair-rule genes (PRGs). One of these genes was thought to be the type-II transmembrane protein encoding gene Tenascin-m (Ten-m (aka odd Oz)). Ten-m, however, does not have a pair-rule function in Drosophila, despite its characteristic PRG-like expression pattern. A recent study in the beetle Tribolium castaneum showed that its Ten-m gene is not expressed like a segmentation gene, and hence is very unlikely to have a function in segmentation.\n\nResultsIn this study, I present data from a range of arthropods covering the arthropod tree of life, and an onychophoran, representing a closely related group of segmented animals. At least one ortholog of Ten-m/odz in each of these species is expressed in the form of transverse segmental stripes in the ectoderm of forming and newly formed segments - a characteristic of genes involved in segmentation.\n\nConclusionsThe new expression data support the idea that Ten-m orthologs after all may be involved in panarthropod segmentation.

Phenotypic plasticity is widely regarded as important for enabling species resilience to environmental change and for species evolution. However, insight into the complex mechanisms by which phenotypic plasticity evolves in nature has been limited by our ability to reconstruct evolutionary histories of plasticity. By using part of the molecular mechanism, we were able to trace the evolution of pre-feeding phenotypic plasticity across the class Echinoidea and identify the origin of plasticity at the base of the regular urchins. The neurosensory foundation for plasticity was ancestral within the echinoids. However, coincident development of the plastic trait and the neurosensory system was not achieved until the regular urchins, likely due to pleiotropic effects and linkages between the two colocalized systems. Plasticity continues to evolve within the urchins with numerous instances of losses associated with loss of sensory capabilities and in one case loss of neurons, consistent with a cost associated with maintaining these capabilities. Thus, evidence was found for the neurosensory system providing opportunities and constraints to the evolution of phenotypic plasticity.

Molluscs are one of the most morphologically diverse clades of metazoans, exhibiting an immense diversification of calcium carbonate structures, such as the shell. Biomineralization of the calcified shell is dependent on shell matrix proteins (SMPs). While SMP diversity is hypothesized to drive molluscan shell diversity, we are just starting to unravel SMP evolutionary history and biology. Here we leveraged two complementary model mollusc systems, Crepidula fornicata and Crepidula atrasolea, to determine the lineage-specificity of 185 Crepidula SMPs. We found that 95% of the adult C. fornicata shell proteome belongs to conserved metazoan and molluscan orthogroups, with molluscan-restricted orthogroups containing half of all SMPs in the shell proteome. The low number of C. fornicata-restricted SMPs contradicts the generally-held notion that an animals biomineralization toolkit is dominated by mostly novel genes. Next, we selected a subset of lineage-restricted SMPs for spatial-temporal analysis using in situ hybridization chain reaction (HCR) during larval stages in C. atrasolea. We found that 12 out of 18 SMPs analyzed are expressed in the shell field. Notably, these genes are present in 5 expression patterns, which define at least three distinct cell populations within the shell field. These results represent the most comprehensive analysis of gastropod SMP evolutionary age and shell field expression patterns to date. Collectively, these data lay the foundation for future work to interrogate the molecular mechanisms and cell fate decisions underlying molluscan mantle specification and diversification.

Arthropods have an incredible diversity of limbs that are modified for walking, chewing, cleaning, mating, grasping, sensing, and more. Understanding the relationships and evolutionary histories of different limbs is a central task, but their sheer diversity makes this a daunting if not impossible task using morphology alone. Here, the in situ expression patterns and CRISPR-Cas9 phenotypes for the five best-studied leg-patterning genes - Distal-less, Sp6-9, dachshund, extradenticle, and homothorax - are described for all limbs of the crustacean Parhyale. Crustaceans are well-suited for this task because their limbs are more diverse than those of other arthropods, and each individual possesses a wide range of limb types that are relevant to many other arthropods, living and extinct. These results will a) provide a template for understanding the genetic basis of limb construction in arthropods more generally based on the strong phenotypes that can be obtained with CRISPR-Cas9, and b) contribute to our understanding of the evolution and affinities of highly modified legs like mouthparts and genitalia using molecular methods to complement previous morphological and embryological approaches.

Asymmetries are essential for proper organization and function of organ systems. Genetic studies in deuterostomes have shown signaling through the Nodal/Smad2 pathway plays a key, conserved role in the establishment of body asymmetries. While Nodal signaling is required for the establishment of left-right asymmetry (LRA) across bilaterian species, little is known about the regulation of Nodal signaling in spiralians. Here, we identified orthologs of the egf-cfc gene, a master regulator of the Nodal pathway in vertebrates, in several invertebrate species, the first evidence of its presence in non-deuterostomes. Our functional experiments indicate that despite being present, egf-cfc does not play a role in the establishment LRA in gastropods. However, experiments in zebrafish suggest that a single amino acid mutation in the egf-cfc gene in the deuterostome common ancestor was the necessary step in inducing a gain-of-function in LRA regulation. This study shows that that the egf-cfc gene likely appeared in the bilaterian stem lineage, before being adopted as a master mechanism to regulate the Nodal pathway and the establishment of LRA in deuterostomes.

Tunicates offer a suitable system to study the evolutionary mechanisms underlying lifestyle transitions, ranging from the biphasic lifestyle of ascidians alternating swimming larva with sessile adults, to the fully free-swimming appendicularians. The appendicularian Oikopleura dioica, despite having only ten pairs of tail muscle cells, roughly half than in ascidian larvae, exhibits a remarkably rich repertoire of behaviours required for house inflation, swimming, nodding, and the rhythmic water pumping through the house. This apparent paradox suggests that functional diversification might have arose not from increased cell number, but from molecular specialization within a minimal musculature. To address this question, here, we reconstruct the evolution of the Myosin class II heavy-chain (Myh) gene family across tunicates and generated a developmental expression atlas of all O. dioica Myh genes. Phylogenetic analysis reveals that the cardio-paraxial Myh subfamily duplicated into two tunicate-specific subfamilies, followed by independent bursts of paralogue duplications in ascidians and appendicularians. In O. dioica, two distinct Myh-Scp-Tb paralogs show cardiac expression, and surprisingly, combinations of different Myh-Scp-Tb paralogues define multiple muscle cell identities along the tail. The innovation of this anteroposterior muscle regionalization associated to the Myh-Scp-Tb expansion provides a plausible mechanism for fine-scale modulation of contractile properties within a minimal musculature, helping to resolve the paradox of "fewer cells, but more contractile properties" in appendicularian tail movements. Conversely, the loss of the Myh-Sj/bw gene, which in ascidians is expressed with post-metamorphic body-wall and siphon muscles, is consistent with a pattern of regressive evolution associated to the absence of trunk muscles in O. dioica. This finding supports the view that the appendicularian lifestyle is a secondary derived condition, and that the last common ancestor of tunicates likely possessed an ascidian-like biphasic lifestyle. Together, our findings indicate that both gene-family duplication and gene loss have shaped the evolution of appendicularian tail muscle, enabling the emergence of complex tail-driven behaviours, and offering new insights into the ancestral lifestyle of tunicates.

The regulation of gene expression is one of the key evolutionary processes driving phenotypic divergence among species. Here, we investigate the tissue-specific gene expression of a non-model adaptive radiation of Cyprinodon pupfishes, characterized by their divergent dietary niches and exceptionally fast rates of craniofacial evolution. By comparing tissue-specific gene expression in the most morphologically divergent skeletal structure, the oral jaws, with the relatively morphologically conserved caudal tail region, we identified genes that were differentially expressed exclusively in the developing jaws of each of the three trophic specialists at hatching (8 dpf) and not in any other species. We then assessed their overlap (as transcriptionally-regulated genes) with adaptive regulatory variants identified in previous genomic studies. Our analysis identified pycr3 and atp8a1 as the most promising for craniofacial evolution in the scale-eaters, both genes with no known previous craniofacial function. We functionally confirmed the craniofacial expression of these genes by in situ mRNA hybridization chain reaction and demonstrated their species-specific expression in branchial and muscle tissues between sister species of this young radiation. Our work underscores the power of integrating tissue-specific transcriptomics with speciation genomics to identify novel craniofacial candidate genes controlling divergent morphogenesis in a natural evolutionary mutant system.

To understand how novel animal patterning emerged, one needs to ask how the development of color patterns has changed among diverging species. Here we examine three species of fruit flies - Drosophila guttifera (D. guttifera), Drosophila palustris (D. palustris), and Drosophila subpalustris (D. subpalustris) - displaying a varying number of abdominal spot rows that were either gained or lost throughout evolutionary time. Through in situ hybridization, we examine the mRNA expression patterns for the pigmentation genes Dopa decarboxylase (Ddc), tan (t), and yellow (y) during pupal development. Our results show that Ddc, t, and y are co-expressed in identical patterns, each foreshadowing the adult abdominal spots in D. guttifera, D. palustris, and D. subpalustris.

Specific character traits of mineralized endoskeletal tissues need to be clearly defined and comprehensively examined among extant chondrichthyans (elasmobranchs, such as sharks and skates, and holocephalans, such as chimaeras) to understand their evolution. For example, tiles of mineralized polygonal structures called tesserae occur at cartilage surfaces in chondrichthyans, but recent studies showing trabecular mineralization at elasmobranch cartilage surfaces suggest that tesserae are not as common as previously thought. Also, while areolar mineralized tissue in elasmobranchs is generally considered a unique, shared chondrichthyan feature, some chondrichthyan species demonstrate bone-like tissues in both a specific region of tesserae termed the cap zone and continuous (not tiled) mineralized neural arches. To clarify the distribution of specific endoskeletal features among extant chondrichthyans, adult skeletal tissues in a holocephalan chimaera (spotted ratfish) and two elasmobranchs (small-spotted catshark and little skate) were characterized using synchrotron radiation and desktop micro-CT imaging, and histological and immunofluorescent assays. Endoskeletal mineralization in the ratfish, catshark, and little skate varied both quantitively in tissue mineral density (TMD), and qualitatively in the morphology and localization of mineralized structures and tissues. For example, TMD of several skeletal elements was significantly lower in ratfish, compared to catshark and little skate. Trabecular and areolar mineralization were shared among these extant chondrichthyan species, but tesserae and bone-like tissues were not. Interestingly, three separate analyses argued that the adult chimaera endoskeleton has features of the embryonic little skate endoskeleton. Generally, this study proposes specific terminology for character states of the extant chondrichthyan endoskeleton and infers those states in ancestral chondrichthyans with reference to fossil data.

Eye evolution has long attracted interest, yet how multiple cerebral eyes within a lineage originate and diversify remains unclear. Annelids display exceptional diversity in eye number and structure, but the homology and function of distinct eye pairs are poorly understood. Here we investigate the cerebral eyes of the sedentary annelid Malacoceros fuliginosus using an integrated developmental, molecular, ultrastructural, and connectomic approach. We show that both the early-developing ventral and later-developing dorsal eyes are simple, few-celled, inverted rhabdomeric eyes that express transcription factors with conserved roles in animal eye development. Two r-opsin paralogs and distinct neurotransmitters are differentially expressed in different photoreceptor cells of the eyes. Ultrastructural reconstructions across larval development reveal differences in cellular composition and growth dynamics, while axonal tracing shows that photoreceptors from ventral and dorsal eyes project to overlapping regions of the larval brain. The overall organization and projections resemble those described in the errant annelid Platynereis dumerilii. Together, these data support the hypothesis that an ancestral cerebral eye duplicated early in annelid evolution, giving rise to multiple eye pairs with stage-specific functions.

The olfactory system plays a critical role in mammalian environmental perception, with some clades relying on an expanded accessory olfactory (vomeronasal) system (VNS) to detect larger, non-volatile odorants. Mammals make extensive use of this system for social communication between conspecifics. Recent studies have begun to investigate how the VNS changes in response to or as part of ecological transitions. Several studies have identified trends of VNS-associated gene loss or regression in secondarily aquatic mammals. However, continuing discussion on genotype-phenotype correlation within the VNS means that greater effort should be made to investigate the morphology of the VNS in species where it remains poorly understood. Here, we use skeletal and soft-tissue data to demonstrate that the vomeronasal groove, an established osteological correlate for the VNO in bats and primates, is also a valid indicator for its presence in Caniformia. Additionally, we confirm the presence of the VNO in the secondarily aquatic North American river otter (Lontra canadensis) and compare its morphology with that of two close-related species, the semi-aquatic American mink (Neogale vison) and the terrestrial long-tailed weasel (Neogale frenata). This study expands the valid taxonomic scope of the vomeronasal grooves proxy as an osteological correlate, confirms the presence of the VNO in the previously undescribed system of the North American river otter, and highlights the complexity of the mammalian accessory olfactory system.

Despite an abundance of gene expression surveys, comparatively little is known about Hox gene function in Chelicerata, with emphasis on the Hox logic of the anterior prosomal segments, which bear the mouthparts. Previous investigations of individual paralogs of labial (lab) and Deformed (Dfd) in the spider Parasteatoda tepidariorum have shown that these play a role in tissue maintenance of the pedipalpal segment (labial-1) and in patterning the first walking leg identity (Deformed-1), respectively. However, broader extrapolations of these data points across chelicerates are hindered by the existence of duplicated copies of Hox genes in arachnopulmonates (e.g., spiders and scorpions), which have resulted from an ancient whole genome duplication event. Here, we investigated the function of single-copy orthologs of lab in the harvestman Phalangium opilio, an exemplar of a lineage that was not subject of this whole genome duplication. Embryonic RNAi against lab resulted in homeotic transformations of pedipalps to chelicerae, as well as reduction and fusion of the pedipalpal segment with adjacent segments. To test for combinatorial function, we performed double knockdown of lab and Dfd, which results in homeotic transformation of both pedipalps and first walking legs into cheliceral identity, whereas the second walking leg is transformed into a pedipalpal identity. Taken together, these results elucidate a model for the Hox logic of head segments in Chelicerata. To substantiate the validity of this model, we additionally performed expression surveys for duplicated copies of lab and Dfd in scorpions and horseshoe crabs, toward understanding the genetic basis of a heteronomous prosoma. We show that repetition of morphologically similar appendages is correlated with uniform expression levels of the Hox genes lab and Dfd, irrespective of the number of gene copies.

Recent molecular phylogenetic studies have raised two questions about the evolutionary history of the calcified exoskeleton of mollusks. The first question concerns the homology of the two types of skeleton; whether spicules and shell plates share an evolutionary origin. The second question is the homology of the shell plates between chitons and other mollusks, including gastropods and bivalves. To gain insight into these questions, we examined the early development of shell plates and spicules in chitons. We identified several developmental genes that are involved in both shell plates and spicules, suggesting that spicules and shell plates share a common evolutionary origin. We also found that subpopulations of the dorsal shell field (the ridge and the plate field) have specific gene expression profiles. The differential gene expression of the ridge and plate field is not identical to the profiles of the zones of the gastropod shell field. This observation may suggest an independent evolutionary origin of the shell plates in chitons and gastropods.